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Abstract. The energy range encompassing the ankle 
of the cosmic ray energy spectrum probably marks 
the exhaustion of the accelerating sources in our 
Galaxy, as well as the end of the Galactic con- 
finement. Furthermore, this is the region where the 
extragalactic flux penetrates the interstellar medium 
and starts, progressively, to be dominant. Although at 
lower energies it is likely that an "average" popula- 
tion of supernova remnants can be defined to account 
for most of the cosmic ray flux, this assumption is 
increasingly difficult to maintain as higher energies 
are considered. One possibility is that supernovas are 
still a main contributor along the first branch of the 
ankle region, but that the acceleration is now coming 
from well localized regions with a characteristic in- 
terstellar medium, or a sub-population of supernovas 
exploding in a peculiar circumstellar environment. 
These possibilities are analyzed in the present work 
using a two-dimensional diffusion model for cosmic 
ray propagation. Special emphasis is given to the 
inner 200 pc of our Galaxy and to the spiral arm 
structure in relation with the Sun position inside the 
disk. 
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I. Introduction 

In our previous work JH, Q, we have analyzed the 
matching conditions of the Galactic and extragalactic 
components of cosmic rays (CR) along the second knee 
and the ankle. From this analysis, it seems clear that an 
acceptable matching of the Galactic and extragalactic 
fluxes can only be achieved if the Galaxy has additional 
accelerators, besides the fiducial Supernova Remnants 
(SNRs) assumed there, operating in the interstellar 
medium. Despite the fact that a different acceleration 
mechanism can be invoked to account for the highest 
energy side of the Galactic spectrum, acceleration by 
SNRs may still play a dominant role. Actually, a possible 
Galactic contribution, dominating at the highest energies, 
could be represented by compact and highly magnetized 
SNRs, like those occurring in the central, high density 
regions of the Galactic bulge, inside the dense cores of 
molecular clouds or those expanding into the circum- 
stellar winds of their progenitors. In the same way, the 
highest energy end of the Galactic spectrum could be 
the result of a non-homogeneous SNR population drawn 
from a spectrum of progenitor masses and evolving in 



different environments corresponding to the various gas 
phases that fill the interstellar medium Q. 
These possibilities are explored in this paper calculating 
with a diffusion Galactic model the contribution of 
localized sources to the Galactic cosmic ray (GCR) 
spectrum, SNRs localized in the Galactic center (GC) 
and in the Galactic Ring presenting the highest CO 
emissivity. These two contributions, along with the 
Galactic spectrum resulting from fiducial SNRs (called 
standard in the following), are combined with a mixed 
extragalactic (EG) spectrum and normalized to match 
HiRes Hi, and Auger 0, Q experimental data. The 
EG model considered is the mixed composition model 
by Allard et al. ID in the case of the uniform source 
distribution model. 

II. SNR EVOLUTION 

The maximum acceleration energy by SNRs is related 
to the shock radius at the end of the adiabatic phase. 
In this section we estimate the radius and age of a 
blastwave at the end of its adiabatic phase for two 
different characteristic cases, a uniform ambient density 
and a power-law density gradient p oc r -2 , the latter 
corresponding to explosions in a pre-existing supersonic 
wind. 

The dynamic of the blastwave in stationary and homo- 
geneous media, is described by the self-similar Sedov- 
Taylor (ST) solution, valid in an ambient with negligible 
pressure Pn — l9l. In media with power-law density 
variations po(r) oc r~ kp (for k p < /c pcr iJ3), the evolution 
of a blastwave of initial energy Eq in the adiabatic phase 
is given by 
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where R s is the radius of the blastwave (normalized to 
some fiducial value J? s (l)), — 1Q ^i° erg and no(l) is 
the ambient hydrogen number density in cm~ 3 at 
The quantity £ in eq. Q] is given by 

t=T^> (2) 
where n — = 2 , and a is the ratio between the total 

' 5-fcp 

energy of the gas inside the blastwave and its kinetic 

l k pcrit = "^rj"> where 7 is the specific heats ratio at the shock. 
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energy. For a specific heat ratio 7 = 
value 

3(5-fc p )(10-3fc p ) 
4 8vr(3 - k p ) 

The age t correspondent to the shock radius R s is 
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,(l-fc p /3)^5i, 
The adiabatic phase ends when radiative cooling dom- 
inates the blastwave evolution. The transition point can 
be estimated as the radius R c (time t c ) at which the half 
of the initial energy has been radiated away: 

l/(7-fc„) 
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where I -1/2 is the radiative form factor for line- 
cooling. The corresponding age of the blastwave t c can 
be estimated using eq. [4] 

From this equation the radius and age of blastwaves at 
the end of the adiabatic phase can be calculated for our 
two cases (for 7 = 5/3): 
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A uniform density ISM is a good approximation for 
SNRs evolving inside the GC region, as well as for the 
main SNR component ("standard") expanding inside 
the general ISM elsewhere in the Galactic disk. In fact 
the ISM inside the inner 200 pc of the Galaxy has a 
density of ~ 10 4 cm 3 , which is very similar to that 
of the the cores of molecular clouds, but extended 
over a large region where several compact remnants 
can be easily observed in radio at different stages of 
evolution iflOl . Furthermore, this region is permeated 
by a magnetic field almost three orders of magnitude 
larger than in the rest of the Galactic disk. A power 
density variation of the form r~ 2 (i.e., k p = 2), on the 
other hand, corresponds to a SN blast wave propagating 
through the circumstellar wind of its progenitor. It is 
unlikely that these circumstellar regions could preserve 
such a well structured density profile inside the high 
density, high turbulence environment of the Galactic 
center. Nevertheless, a k p — 2 profile could described 
well the ambient medium of SNR expanding inside the 
still dense, but rather isolated regions, associated with 
dispersed molecular clouds inside the Ring. 

The maximum acceleration energy for a particle of 
charge Z is proportional to ZeBR, where B and R 



are the magnetic field intensity and the size of the 
acceleration region. If we consider a standard ISM and 
the GC, typical values for gas density and magnetic 
field are no ~ 1 cm~ 3 and B ~ pG, for the first and 
no ~ 10 4 cm -3 and B ~ mG for the second. The 
maximum scale R associated with particle acceleration 
can be estimated as the SNR radius at the transition point 
between the adiabatic and the radiative phase, given by 
eq. [6] or [8] depending on the value of k p . 
The ratio between the maximum energy achievable by 
protons by acceleration in SNRs in the two cases is: 
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In the case of uniform gas density, using eq. [6] the ratio 
becomes 
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III. Diffusion Galactic model 



We used the numerical diffusive propagation code 
GALPROP lfT2l . |[T3l to reproduce the galactic spectrum 
from SuperNova Remnants (SNRs). The diffusive model 
is axisymmetric. The propagation region is, in cylindrical 
coordinates, bounded by R = Rh = 30 kpc and 
z = Zh = 4 kpc, beyond which free escape is assumed, 
where R is the Galactocentric distance and z is the 
altidude from the Galactic plane. 
The propagation equation is: 



dtp 



-£ = q(r,p) + W '(AcsVVO 
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where ijj{f,p,t) is the density per unit of total particle 
momentum, q(f,p) is the source term, D xx is the spatial 
diffusion coefficient, p = dp/dt is the momentum loss 
rate and 77 and r r are the time scale of fragmentation 
and the time scale of radioactive decay respectively. The 
diffusion coefficient is taken as (3Do(p/ pd) 6 , where p 
is the particle rigidity, Do is the diffusion coefficient at 
a reference rigidity po and 5 = 0.6. 
Stable nuclei with Z < 26 are injected, at the source 
with energy independent isotopic abundances derived 
from low energy CR measurements fl3l . 
Detailed and realistic interstellar molecular (Ha), atomic 
(H) and ionized (HI) hydrogen distributions are used 

m. 

The distribution of cosmic rays sources used for stan- 
dard SNRs is that of Galactic SNRs deduced from 
EGRET gamma-ray data [12], while for SNRs operating 
in the GC and in the Ring we used a uniform distribution 
of sources limited to the regions R < 0.2 kpc, \z\ < 0.2 
kpc and 4.5 < R < 5.5 kpc, \z\ < 0.1 kpc, respectively. 
The injection spectrum is a a power law function in 
rigidity with a break at rigidity po, beyond which it falls 
exponentially with a rigidity scale p c : 
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I{P) = 




exp 
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where a = 2.05, while po and p c depend on the SNRs 
considered: 

. standard SNRs: pg* = 3.6 PV, pf = 2.5 PV; 

. Ring SNRs: p R = 17.5 PV, pf = 12.15 PV; 

. GC SNRs: p^ c = 80 PV, p^ c = 55.6 PV. 

The rigidity scale and rigidity cut-off for standard 
SNRs were determined empirically so that the total GCR 
spectrum fits the position of the first knee measured by 
KASCADE (E knee ~ 5 PeV) [ 1 1 1 and the shape of the 
observed spectrum beyond it up to the highest possible 
energies. The estimate of the rigidity cut-off for the 
SNRs in the GC and in the Ring were also determined 
empirically but in order to fit the measured total spec- 
trum, once combined with the standard SNRs spectrum 
and with the EG spectrum. This estimate was performed 
taking into account the estimate of E~~ iX / E^ ax given 
by eq. HU which suggests a rigidity cut-off for the GC 
component Pq C ~ 20 x p$. 

The physical conditions required for the Ring region in 
order to provide the rigidity cut-off p R = 17.5 PV can 
be estimated from eq. [8] 

The ratio between the maximum energy achievable by 
protons in the case of SNRs exploding in a pre-existing 
wind in the Ring region and standard SNRs is: 



ipRing 
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For a magnetic field density <~ lQpG characteristic of 
the Ring region and the typical values of the Galactic 
ISM, using eqs.|6][8]and our empirical ratio p R / pg* ~ 5, 
the required gas density in the Ring region is of order 
of 200 cm~ 3 . This density value is reasonable for the 
Ring region, characterized by heterogenous medium of 
regular ISM interspersed with molecular cloud cores of 
density as high as 10 4 -10 5 cm~ 3 . 

IV. Diffusive Galactic spectrum 

The calculated diffusive Galactic spectra from the 
three sources are combined with the EG spectrum and 
normalized in order to match KASCADE data at E ~ 
3 x 10 6 GeV fl~Q and HiRes data 0, at higher 
energy. The resulting CR spectrum is shown in Fig. Q] 
superimposed to several experimental data results. 
We calculated the number of sources contributing to 
the different Galactic flux components. Integrating the 
source distribution functions and assuming that the total 
energy pumped by a SNR into the CR component is a 
constant fraction of its total kinetic energy independently 
on its environment, we found that ~ 6500 standard 
SNRs are required in order to account for the main 
component of the CR Galactic flux up to the second knee 
energy, while the fluxes coming from the GC and from 
the Ring requires ~ 7 and ~ 260 sources, respectively. 



The number of required SNRs in the GC region is 
sufficiently small to be readily supplied by the observed 
population inside the inner 200 pc of our Galaxy. In fact, 
such a small number could pose a potential problem, 
since random variations in the number of SNRs could 
distort the shape of the spectrum along the ankle as a 
function of time. 

In order to check the extent of their possible effect, we 
have perturbated the number N of SNRs in the GC 
population by ±\/~N (and analogously for the number 
of sources in the Ring). The results are shown in Fig. 
[T] by the curves limiting the hatched areas, where it 
can be seen that the ankle remain smooth but, under 
extreme conditions, its location could change in energy 
by up to half a decade under stochastic fluctuations of 
the population of high energy SNRs. In any case, since 
the average time between SN in the GC is much smaller 
than the 10 4 yr of diffusion time up to the solar circle, 
the time dependence of the ankle should be further 
suppressed. 

The same procedure of normalization of the Galactic 
components has been applied in order to match the total 
spectrum with Auger data 0, at high energy and to 
KASCADE data at E ~ 3 x 10 6 GeV ifTTl . 
The resulting total spectrum and each Galactic contribu- 
tion are shown in Fig. [2] In this case, the CR flux from 
sources in the GC is considerably lower with respect to 
the HiRes case, while the Ring and standard components 
are of the same order of magnitude. Integrating the 
source distribution functions, the number of sources 
contributing to the GCR flux is ~ 6900 for standard 
SNRs and ~ 230 and ~ 4 for the Ring and GC, 
respectively. As in the previous case, the result for a 
perturbation of the number of sources in the GC region 
and in the Ring are indicated by the curves limiting the 
hatched area [2] Even if the number of sources in the GC 
are smaller with respect to the HiRes case, the global 
behavior is consistent with that observed in the HiRes 
case. 

V. Conclusion 

A previous analysis 0, Q suggests that additional 
acceleration mechanisms are required besides accelera- 
tion from fiducial (standard) SNRs in order to account 
for GCR flux up to the highest energy. In a scenario 
where accelerators different from SNRs are disregarded, 
the highest energy part of the galactic spectrum could 
be the result of SNRs evolving in peculiar environments. 
We explore here the possible contribution of SNRs inside 
the GC and in the dense Galactic Ring localized at 
4.4 < R < 5.5 kpc. We have estimated empirically the 
maximum energy, rigidity cut-off and rigidity scale re- 
quired of these additional components in order to fit the 
total measured spectrum. Using the Galactic diffusion 
model GALPROP, we have verified that these parameters 
are compatible with the evolutionary properties and 
expected CR luminosities of SNR in these two regions. 
Therefore, the present analysis suggests that acceleration 
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Fig. 1. Diffusive total Galactic spectrum combined with the mixed EG spectrum (&eg) ar| d normalized to KASCADE [111 and HiRes 

data (4), 0. The different contributes from standard SNRs (<I>g), SNRs in the GC (&q C ) and in the Ring (<I>g) are shown. The curves 
limiting the hatched areas correspond to a variation ±\/~N of the number of contributing sources in the GC and in the Ring. HiRes data and 
other several experimental data results are shown. 




Fig. 2. Diffusive total Galactic spectrum (<E>g*) combined with the mixed EG spectrum ($bg) an d normalized to KASCADE [11] and Auger 
data. The different contributes from standard SNRs SNRs in the GC ($g c ) and in the Ring ($g) are shown. The curves 

limiting the hatched areas correspond to a variation ±V~N of the number of contributing sources in the GC and in the Ring. Auger data and 
other several experimental data results are shown. 



by SNRs has the potential to account for the whole 
Galactic cosmic ray flux if three general populations 
are considered: (i) a main population of standard SNRs 
which evolve in media with no ~ 1 cm~ 3 and B ~ 1/iG, 
(ii) - 2 x 10 2 SNRs associated with ISM of n Q ~ 10 2 
cm~ 3 and B ~ inside the inner Ring and (iii) 

< 10 SNRs immersed in the GC ISM, where n ~ 10 4 
cm~ 3 and B ~ lmG 
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